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ABSTRACT 

Context. We report on near-infrared (IR) observations of the three anomalous X-ray pulsars XTEJ18 10-197 , IRXS 11708-^009 and 
IE 1841-045 , and the soft gamma-ray repeater SGR 1900+14 , taken with the ESO-VLT, the Gemini, and the CFHT telescopes. 
Aims. This work is aimed at identifying and/or confirming the IR counterparts of these magnetars, as well as at measuring their 
possible IR variability. 

Methods. In order to perform photometry of objects as faint as K,, ~ 20, we have used data taken with the largest telescopes, equipped 
with the most advanced IR detectors and in most of the cases with Adaptive Optics devices. The latter are critical to achieve the sharp 
spatial accuracy required to pinpoint faint objects in crowded fields. 

Results. We confirm with high confidence the identification of the IR counterpart to XTEI18 10-197 , and its IR variability. For 
IE 1841-045 and SGR 1900+14 we propose two candidate IR counterparts based on the detection of IR variability. For IRXS J 1708- 
4009 we show that none of the potential counterparts within the source X-ray error circle can be yet convincingly associated with 
this AXP 

Conclusions. The IR variability of the AXP XTEJl 8 10-197 does not follow the same monotonic decrease of its post-outburst X-ray 
emission. Instead, the IR variability appears more similar to the one observed in radio band, although simultaneous IR and radio 
observations are crucial to draw any conclusion in this respect. For IE 1841-045 and SGR 1900+14 follow-up observations are 
needed to confirm our proposed candidates with higher confidence. 

Key words. stars:pulsars:individual:XTEI1810-197 , IRXS J 1708-4009 , IE 1841-045 , SGR 1900+14 



1. Introduction 

■ In recent years, the study of isolated neutron stars has become 
one of the most challenging research areas in high-energy astro- 
physics, largely as a result of the discovery of several new classes 
of sources besides the classical, well-studied radio pulsars. The 
most extreme objects are the so-called "magnetars". This class 
comprises the Anomalous X-ray Pulsars (AXPs) and the Soft 
Gamma-ray Repeaters (SGRs), observationally very similar in 
many respects. They all are slow X-ray pulsars with periods in 
a narrow range (P= 2-12 s), relatively large period derivatives 
(P = lO-i"* - lO^'^ss-i), spin-down ages of 10^^ - lO^'yr, and 
magnetic fields, as infeiTed from the classical magnetic dipole 
spin-down formul^ of lO''* - 10'^ G, larger than the electron 
quantum critical field - 4.4 x 10'^ G) above which quan- 
tum effects become crucial. AXPs and SGRs are strong X-ray 
emitters, with X-ray luminosities of about 10-'"' - lO-'^ergs"'. 
Their 0.1-lOkeV persistent emission has relatively soft spectra 
usually modeled by an absorbed blackbody (kT~ 0.2-0.6) plus a 
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power-law (F ~ 2^; for a review see e.g. Woods & Thompson 
2006). Their X-ray energy output is much larger than the rota- 
tional energy loss and implies that these sources are not rotation- 
ally powered, at variance with most young isolated neutron stars. 
Rather, AXPs and SGRs are believed to be powered by the neu- 
tron star ultra-strong magnetic field (Duncan & Thompson 1992; 
Thompson & Duncan 1995). In the "magnetar" model, crustal 
deformations, driven by magnetic stresses imparted to the crust 
by the strong internal toroidal magnetic field, are responsible for 
the observed activity, X/y-ray bursts and giant flares. 

Until not long ago, the persistent emission of AXPs and 
SGRs was detected only in the soft X-ray range. However, in 
the last few years the availability of more sensitive optical and 
IR telescopes, as well as of y-ray satellites, opened new windows 
in the study of AXPs and SGRs, unveiling their multi-band prop- 
erties. In particular, observations from ESO, Gemini, CFHT and 
Keck (see Israel et al. 2004a for a review) led to the discovery 
of faint (K5 ~20), and in some cases variable, IR counterparts 
to five out of seven confirmed AXPs (one of which was also 
detected in the optical band by HuUeman et al. (2000)), and to 
one (out of four) confirmed SGR jlsrael et al.l2005LlKosugi et al.l 
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Telescope/Instrument 


Target 


Date (YY.MM.DD) 


MJD 


Filter 


Exposure Time (s) 


FWHMC) 


airmass 


VLT/NACO 


XTEJ181()-197 


2003.10.08 


52920 


H 


2520 


0.189 


1.17 






2003.10.08 


52920 


Ks 


2520 


0.243 


1.49 






2004.03.14 


53078 


H 


1080 


0.108 


1.47 






2004.03.14 


53078 


K,s 


1200 


0.094 


1.20 






2004.03.14 


53078 


J 


1080 


0.081 


1.13 






2004.09.09 


53257 


J 


1800 


0.170 


1.01 






2004.09.09 


53257 


H 


3600 


0.140 


1.08 






2004.09.13 


53261 


Ks 


3600 


0.099 


1.01 




IRXS J 1708^009 


2003.05.20 


52779 


K,s 


2400 


0.083 


1.04 






2003.06.19 


52809 


H 


2400 


0.122 


1.40 






2003.06.20 


52810 


J 


2400 


0.122 


1.04 






2004.03.25 


53089 


L' 


2100 


0.119 


1.04 




SGR 1900+14 


2006.04.01 


53826 


K,s 


1560 


0.125 


1.43 






2006.07.21 


53937 


Ks 


2040 


0.111 


1.21 


Gemini/NIRI 


XTEJ1810-197 


2003.09.18 


52900 


Ks 


285 


0.397 


1.62 






2003.09.18 


52900 


H 


285 


0.397 


1.99 






2003.09.18 


52900 


J 


285 


0.397 


1.78 






2004.06.10 


53166 


Ks 


510 


0.631 


1.32 






2004.06.10 


53166 


H 


510 


0.526 


1.30 






2004.06.10 


53166 


J 


510 


0.655 


1.32 






2004.07.28 


53214 


Ks 


510 


0.468 


1.39 






2004.07.28 


53214 


H 


510 


0.526 


1.53 


CFHT/AOBIR 


IE 1841-045 


2002.08.17 


52503 


K' 


1260 


0.137 


1.11 






2002.08.17 


52503 


H 


1260 


0.161 


1.31 



Table 1. Observations summary. 




Fig. 1. XTE J1810-197 : Left panel: VLT/NACO K-band image (March 2004). The radio position jCamilo et al.ll2007cl) is marked 
together with its candidate counterpart (object A). Twelve nearby objects, some of which used as a reference for relative photometry, 
are also marked in the figure. Right panel: ¥^s vs.H-Ks colour-magnitude diagram of the sources detected in the NACO field of view 
(filled circles). Object A is marked with a rectangle indicating the range in colour and magnitude covered throughout the epochs, 
while some of the nearby objects marked in the map and listed in Tab.|2]are plotted in red and labeled. Magnitudes and colours for 
all the sources but object 'A' have been averaged over all the observations. 



l2005h . Furthermore, deep INTEGRAL observations in the 20- 
200 keV band revealed that most of these highly magnetize d 
sources are also hard X-ray emitters dKuiper et al.112004 l2006h . 
The recent discovery of pulse d radio emission from the transient 
magnetars XTEJ1810-197 dHalpern et all I2005L ICamilo et alj 
I2OO6 ) and IE 1547.0 -5408 (Camilo et al. 2007a), the only cases 
so far (Burgay et al,Jl2006, 2007), is intriguing since radio emis- 
sion was believed to be quenched in magnetic fields above the 
quantum critical limit (Baring & Harding 1998). 

Despite the more complete observational picture we have 
now for AXPs and SGRs, the physical processes at the basis of 
their emission in the different bands are not fully understood as 
yet. This is particularly true as far as the optical and IR emission 



of these sources is concerned. The extrapolation of the canonical 
blackbody plus power-law, used to model the soft X-ray emis- 
sion of these sources, largely over-predicts their optical and IR 
emission. Recently, new attempts have been made at modeling 
the multi-ban d spectra o f magnetars to overcome this problem. 
In particular, Rea et al.l ( 1200 7 a) showed that AXP spectra can 
be fitted by a resonant cyclotron scattering model or two log- 
parabolic functions, which do not over-predict their dim optical 
and IR magnitudes. However, despite solving this issue, these 
new spectral models are still far from giving any physical inter- 
pretation for the optical and IR emission. 

Some possibilities are that either the magnetars opti- 
cal/IR emission arise from the star surface/magnetosphere 
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(Beloborodov & Thompson 2006), or from the reprocessing of 
the X-ray emission via a fossil disk around the neutron star 
(Chatterjee et al. 2000, Perna et al. 2000). I ndeed, the Spitzer 
mid-IR detection of the AXP 4U0142+614 dWang et al.ll2"006l) 
has been interpreted as the first evidence of a passive fossil disk 
around an AXP, believed to emit through reprocessing of the X- 
ray radiation from the magnetar 

In this paper we report on Gemini and ESO-VLT observa- 
tions of the transient AXPXTE J1810-197 , on ESO-VLT obser- 
vations of the AXP IRXS 11708^009 and of SGR 1900+14, 
and CFHT observations of the AXP IE 1841-045 (§|2ll. The re- 
sults are presented and discussed in the context of different emis- 
sion scenarios, and compared with the optical and IR properties 
of other magnetars (see §|4]i. 

2. Observations and data reduction 

2.1. Observations 

Observations were performed with the 8 m Very Large Telescope 
(VLT) and the 3.5m Canada-France-Hawaii Telescope (CFHT), 
both equipped with adaptive optics (AO) and near-IR cameras. In 
addition, we used archival observations performed with the 8 m 
Gemini telescope, also equipped with an IR camera but with- 
out AO, that were offered starting from 2004 only. A summary 
of the telescopes and instruments used, as well as the observa- 
tions logs, is reported in Tab.[T] In all cases images were obtained 
through the standard observing technique commonly applied to 
IR aiTays, i.e. by stacking repeated exposures along each node 
of a pre-defined dithering pattern, with each node jittered by a 
few arcsec with respect to the previous one. This allows to ob- 
tain a set of dithered exposures that are used to produce sky- 
background images. 

VLT/Yepun observations of XTE J1810-197, IRXS J1708- 
4009 and SGR 1900+14 were carried out in different observing 
runs (see Tab. [T) using NAos COnica (NACOfI with the S27 
camera (0'.'027 pixel size, 28" x 28" field of view). Gemini- 
North observations of XTEJ1810-197 were carried out using 
the Near-InfraRed Imager (NIRlfl with the f/6 camera (0'.'117 
pixel scale, 120" x 120" field of view). IE 1841-045 was ob- 
served at the CFHT using the Adaptive Optics Bonnet (AOBj3 
InfraRed camera (0'.'035 pixel size, 36" x 36" field of view). 
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Fig. 2. XTE Jl 8 10-197 : Top panel: X-ray flux (shown in units of 
10"" ergcm"^ s"') history in the 0.5-lOkeV energy range (from 
Gotthelf & Halpern 2007). Lower panels: from top to bottom: 
Ks, H, and J band relative magnitude variation for XTE J1810- 
197 counterpart (red star) and for two nearby, bright compari- 
son objects. Magnitude variations are computed relative to the 
September 2004 observation. To better fit in the plot, magnitude 
differences for objects #2 and #8 have been shifted by + 1 . 1 and - 
1 . 1 mag, respectively. Bottom panel shows the H-K, colour vari- 
ation. Magnitude eiTors at (Icr) are shown. These include inter- 
nal errors, aperture correction errors and zero point errors. 



2.2. Data reduction and calibration 

The VLT data were pre-reduced with the ESO NACO pipelin^E 
which is based on the package eclipse, while the native version 
of the eclipsJ^puck^age was used for the Gemini and CFHT data. 
In all cases, the image pre-reduction produced a final co-added, 
sky-subtracted and flat-fielded image for each band. 

Depending on the case, fluxes were computed either through 
PSF or aperture photometry ti sing the IRAlQ version of the 
daophot package (IStetsoniri992l) . For AO observations we found 
some instrument-dependent effects which hampered the PSF fit- 
ting and discouraged the use of PSF photometry. In particular, 

2 www.eso.org/instruments/naco 

3 www.gemini.edu/sciops/instruments/niri 

4 www.cfht.hawaii.edu/Instruments 

5 www.eso.org/instruments/naco/ 
* www.eso.org/eclipse 

^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



in the VLT-NACO images the PSF was found to be variable 
across the field of view, while in the CFHT-AOB images it was 
found to feature an unusual profile with a clearly visible sec- 
ondary diffraction ring. Both effects are difficult to account for, 
even by using more sophisticated and position-dependent PSF 
models. On the other hand, in the Gemini-NIRI images the PSF 
was found to be more stable, although worse sampled due to 
the larger pixel size. For these reasons, we decided as a general 
strategy to use aperture photometry for our AO images, i.e. the 
NACO and AOB ones, and to use PSF photometry for the NIRI 
images. For the former, aperture photometry was performed us- 
ing an aperture or diameter twice the measured PSF. Aperture 
and PSF photometries have been compared by applying aperture 
coiTections using the IRAF daogrow algorithm, which calculates 
the aperture correction extrapolating to an ideally infinite aper- 
ture. 

For the VLT images the photometric calibra tion was per- 
formed by observing standard stars from the 'Persso n et al.l 
(1998) catalog, usually observed under similar conditions (and 
at about the same epoch) as our targets. For all the NACO runs 
the zeropoint uncertainty has been found to be very small, i.e. in 
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Source Id. 


R.A.(2000) 


Dec. (2000) 


52900 


52920 


53078 


53166 


53214 


53257 


53261 






272.46283 


-19.731 10 


20.81±0.08 


20.82±0.09 


2 1 . 1 8± 0.05 




21.21±0.I0 




21.04± 0.09 




272 46390 


-19.73016 


13.87±0.07 


13.85±0.04 


14.26± 0.04 




14.01 ±0.04 




14.01± 0.08 


2 




-19.73061 


1 7 23+0 07 


17 2^+0 04 


1 7 23+ 04 




1 7 36+0 04 




1 7 ^0+ 08 


3 


272.46338 


-19.73070 


19.09±0.07 


19.21±0.05 


19.47± 0.04 




19.15±0.05 




19.41± 0.08 


4 


272.46335 


-19.73110 


19.43±0.07 


19.37±0.05 


19.46± 0.04 




19.42±0.04 




19.47± 0.08 


5 


272.46399 


-19.73098 


19.52±0.07 


19.77±0.06 


19.67± 0.04 




19.62±0.04 




19.70± 0.08 


6 


272.46384 


-19.73145 


18.61 ±0.07 


18.61±0.04 


18.60± 0.04 




18.66±0.04 




18.67± 0.08 


7 


272.46327 


-19.72952 


17.69±0.07 


17.70±0.04 


17.73± 0.04 




17.75±0.04 




17.75± 0.08 


8 


272.46535 


-19.72897 


15.48±0.07 


15.40±0.04 


15.45± 0.04 




15.51±0.04 




15.53± 0.08 


9 


272.46264 


-19.72882 


15.25±0.07 


15.15±0.04 


15.23± 0.04 




15.27±0.04 




15.23± 0.08 


10 


272.46082 


-19.73013 


14.43±0.07 


14.19±0.04 


14.55± 0.04 




14.44±0.04 




14.29± 0.08 


11 


272.46430 


-19.73252 


14.83±0.07 


14.74±0.04 


14.91± 0.04 




14.87±0.04 




14.90± 0.08 


12 


272.46502 


-19.73141 


16.59±0.07 


16.50±0.04 


16.51± 0.04 




16.63±0.04 




16.54±0.08 


H 


A 






21.67±0.12 


21.81±0.16 


22.55± 0.07 


21.48± 0.12 


21.73± 0.11 


22.29±0.18 




1 






14.27±0.05 


14.20±0.04 


14.39± 0.04 


14.24± 0.09 


14.24± 0.08 


14.24±0.08 




2 






17.66±0.05 


17.68±0.04 


I7.66±0.04 


17.69± 0.09 


17.71±0.08 


17.67±0.08 




3 






20.00±0.06 


20. 11 ±0.05 


20.07± 0.04 


19.90± 0.09 


19.78± 0.09 


20.07±0.08 




4 






22.24±0.14 


22.05±0.11 


22.06± 0.09 


21.66± 0.23 


21.87± 0.10 


22.64±0.17 




5 






22.28±0.28 


22.74±0.42 




21.79± 0.21 


2I.68± 0.10 


22.80±0.11 




6 






21.17±0.09 


21.38±0.07 


21.38±0.05 


2I.I1±0.10 


2I.15± 0.09 


21.34±0.09 




7 






18.16±0.06 


18.19±0.04 


18.14± 0.04 


18.14± 0.09 


18.15± 0.08 


18.20±0.08 




8 






18.95±0.06 


I9.00±0.05 


I8.90± 0.04 


18.90± 0.09 


18.93± 0.08 


18.95±0.08 




9 






18.42±0.06 


18.42±0.04 


18.34± 0.04 


18.44± 0.09 


18.41±0.08 


18.44±0.08 




10 






14.62±0.05 


14.5 1 ±0.04 


14.65± 0.05 


14.61±0.09 


1 4.61 ±0.08 






11 






15.15±0.05 


15.12±0.04 


15.19± 0.04 


15.13± 0.09 


15.14± 0.08 


15.12±0.08 




12 






16.91±0.05 


16.89±0.04 


16.83± 0.04 


16.93± 0.09 


16.91±0.08 


16.90±0.08 




J 


A 






22.92± 0.22 




23.89 ± 0.23 






23.45 ±0.12 




1 










15.37 ±0.06 






15.04 ± 0.04 




2 






18.85±0.05 




18.83 ±0.06 






18.88 ± 0.04 




3 






21.06± 0.06 




21.35 ±0.07 






21.35 ±0.04 




7 






19.34± 0.05 




19.39 ±0.06 






19.36 ± 0.04 




9 






18.06± 0.04 




18.00 ±0.06 






18.01 ± 0.04 




10 










15.54 ±0.06 










11 






15.98± 0.04 




16.01 ±0.06 






15.97 ± 0.04 




12 






I9.04±0.05 




18.94 ±0.06 






18.99 ±0.04 





Table 2. XTE J1810-197 : From top to bottom: K,, H and J magnitudes for the candidate IR counterpart of XTEJ1810-197 and 
some of the twelve nearby objects, for the 6 different observing runs (see TabO. Numbers atop the magnitude columns give the 
observations epochs in Modified Julian Days (MJD). Magnitude errors are given at Icr and include zero point errors, either absolute 
or relative, that are, for the brightest objects, the dominant error source. 



the 0.02 - 0.05 magnitude range. The residual uncertainty intro- 
duced by the extinction term is almost wiped out by the closeness 
in time and airmass of target fields and standards. Since it turns 
out to be $ 0.01 magnitudes it has been neglected. 

The Gemini images of XTE J18 10-197 have been at first 
instance calibrated with the available standard stars, but the un- 
certainty in the zero points and the zero point differences usu- 
ally found between standards taken on the same nights sug- 
gested us to discard the primary standard calibration. A photo- 
metric calibratio n was also tried using stars from the 2MASS 
jSkrutskie et alj [2006) catalog identified in the target fields. 
However, such an on-the-frame calibration suffered from two 
main problems. Firstly, the number of suitable 2MASS stars is 
often quite small and the brighter ones were found to be gen- 
erally saturated in our images. Furthermore, the 2MASS catalog 
was built from images with a 2" pixel size so that, in many cases, 
the measured magnitudes correspond to blended objects which 
are, instead, resolved in our higher resolution images. Thus, 
given the accuracy of the NACO zeropoints and the very low 
r.m.s. of the photometric solution, we chose to re-calibrate the 
Gemini photometry of the XTE J18 10-197 field on the NACO 
one by using, as secondary photometric calibrators, a set of well- 
suited field stars in common between the NACO and the NIRI 
images. As a reference, we chose the September 2004 NACO 
observations because of their overall better quality. In order to 
ensure full consistency among observations performed at differ- 
ent epochs (Tab. [B, also the NACO October 2003 and March 
2004 photometry of the XTE J1810-197 field was re-calibrated 



to the September 2004 one. The relative photometric uncertainty, 
based on about 100 stars in common among the various images 
is of the order of ^ 0.05 magnitudes for the NACO-to-NACO 
registration, and slightly worse for the Gemini-to-NACO one (» 
0.08 magnitudes). 

For the CFHT photometry, for which no NACO images are 
available, the only option was to use the 2MASS catalog, after 
carefully accounting for all the caveats mentioned above. For 
the L' observation of IRXS J 1708-4009 we have used standard 
stars from the UKIRT IR standard catalog (Leggett et al. 2003|), 
yielding a photometric accuracy of ~ 0.1 magnitudes. 

In all cases, astrometric calibration was performed using 
2MASS stars as reference, yielding an average uncertainty of 
~ 0'.'2, after accounting for the intrinsic astrometric accuracy 
of 2MASS (~ 200 mas) and the r.m.s. of the astrometric fits. 
However, due to the small instrument pixel scales, the latter has 
been found to be much smaller than the 2MASS astrometric ac- 
curacy and hence it has been neglected. 



3. Data analysis and results 

3.1. XTE J1 810-197 

Fig.[I](left) shows the XTE J1810-197 field (K, band) observed 
with NACO on October 2003 (see also Israel et al. 2004b). The 
XTEJ1810-197 candidate IR counterpart, labeled A, is coin- 
cident with the radio position (RA 18'' 09" 5 P. 087, DEC - 
19° 43' 51'.'93; error = tt'Ol; Camilo et al. 2006). For both 
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the NACO and NIRI observations, photometry was computed 
and caUbrated as described in the previous section. For each in- 
strument, the resuhing J, H and band catalogues were then 
matched and compared with the muhi-band ones. The muhi- 
band magnitudes of the XTE J1810-197 candidate counterpart, 
all calibrated to the NACO system, are reported in Tab.|2] to- 
gether with those of a number of field objects, chosen among 
those used as a reference for the photometry cross-calibrations 
(see previous section). We note that the candidate counterpart is 
not detected in the J and K^-band Gemini observations of the 
June 2004 run which were performed in non-photometric con- 
ditions. Given the short integration time, the overall poor data 
quality, and the lack of standard observations, we could only de- 
rive tentative estimates on the limiting magnitudes (3cr upper 
limits of J>21.9 and >20) which are of no use to constrain 
the flux evolution of our source. 

The time-averaged (K,, H-Kj) colour-magnitude diagram 
(CMD) of the field is shown in Fig.[T] (right). A well-defined 
main sequence group of stars, probably tracking a young star 
cluster, is recognizable in the left part of the diagram, while 
a second red clump of stars, most likely red giants or highly 
absorbed main sequence stars, is present on the right. The 
XTEJ1810-197 counterpart (green star) appears much fainter 
than the nearby comparison stars (red dots), and its colour in- 
deed suggests that it might be a peculiar object. 

The variabiHty of the XTEJ1810-197 IR counterpart is 
shown in Fig. ID together with the X-ray variability of the X- 
ray source (see ^for details). While in the X-rays the source 
features a clear, almost monotonic flux decay, in the IR there 
is no clear trend in the flux evolution, and hence no correla- 
tion with the X-rays. Indeed, the flux seems to vary in an er- 
ratic way in all the three IR bands. In particular, we found that 
the overall variability is apparently larger in the H and J bands 
(1.07+0.24 and 0.97±0.32 magnitudes, respectively) than in the 
Kj one (0.40±0.13 magnitudes). We note that our detailed re- 
analysis of the October 2003 and March 2004 observations of 
XTEJ1810-I97 yields a band variability slightly smaller 
(0.36+0.10) than that originally reported by Rea et al. (2004) 
using the same observations (0.50+0.10), although it is still con- 
sistent within the errors. We attribute this (non significant) dif- 
ference to the more accurate relative photometry between the 
two epochs, performed using a larger set of secondary photo- 
metric calibrators than the one used in Rea et al. (2004). The 
last panel shows the (H-K,) colour variability of the source. In 
general, we note that the source seems to become redder when 
its flux decreases, although no clear trend can be recognized in 
our sparse photometry. We need to warn, however, that any con- 
clusion on the source flux and colour variability must be taken 
with due care. First of all, the apparent lack of correlation be- 
tween the variability in the H and bands is only suggested 
by the July 2004 Gemini band observation. If we take this 
observation out and we consider only the co-eval flux measure- 
ments, the K5 band lightcurve becomes more consistent with the 
H (as well as with the J) band one. More generally, although spe- 
cial care was devoted to check the cross-instrument calibration 
(see §2.2), it is possible that our relative photometry is affected 
by random errors. These may be induced, e.g. by fluctuations in 
the atmospheric conditions, sky background, seeing, quality of 
the primary calibration frames (dark and flatfields), and glitches 
in the detector performance. Of course, while these effects are 
marginal for the relatively bright stars that we have used as a 
reference to compute our cross-instrument photometric calibra- 
tion, they are indeed significant for much fainter stars, such as 
the XTEJ1810-197 candidate counterpart (K, ~ 21), and can 



Id 


RA(°) 


DEC (°) 


J 


H 


Ks 


1 


257.19551 


-40.14796 


20.93±0.09 


18.60±0.06 


17.53±0.04 


2 


257.19510 


-40.14815 


22.23±0.17 


19.84±0.09 


18.81±0.05 


3 


257.19560 


-40.14799 


22. 11 ±0.25 


20.01±0.12 


18.85±0.05 


4 


257.19503 


-40.14792 


22.03±0.17 


20.19±0.09 


19.63±0.07 


5 


257.19498 


-40.14800 




20.50±0.10 


19.95±0.08 


6 


257.19495 


-40.14800 






20.03±0.09 


7 


257.19494 


-40.14793 




20.90±0. 1 1 


20.25±0.08 


8 


257.19563 


-40.14790 






20.83±0.10 


9 


257.19502 


-40.14786 




21.55±0.18 


20.92±0.11 


10 


257.19531 


-40.14795 






20.94±0.11 


11 


257.19500 


-40.14804 






20.69±0.10 


12 


257.19548 


-40.14812 






21.12±0.19 


13 


257.19526 


-40.14761 


- 


- 


21.23±0.14 


14 


257.19545 


-40.14774 






21.34±0.21 


15 


257.19523 


-40.14807 






21.52±0.17 


16 


257.19509 


-40.14787 






21.58±0.12 


17 


257.19510 


-40.14799 






21.59±0.16 


18 


257.19538 


-40.14816 




21.34±0.15 


21.66±0.14 


19 


257.19550 


-40.14813 






21.74±0.16 


20 


257.19530 


-40.14792 






21.88±0.19 



Table 3. IRXS J 1708-4009 : J, H, and magnitudes for all the 
objects detected within 1'.' 1 of the source X-ray position, ordered 
by decreasing Kj magnitude. No object has been detected in the 
L' band. Magnitude errors are given at Icr, and include zero point 
uncertainties. 




Fig. 3. IRXS 11708^009 : NACO Kj-band image of the source 
field with the I'.'l radius (99% confidence level) X-ray error cir- 
cle overlaied. Objects detected at > 3cr are labeled. 



increase the overall uncertainties on the measured magnitudes. 
To quantify these random errors we have cross-checked the pho- 
tometry of a number of test stars with brightness comparable to 
the one of our target. We found that all test stars show a scatter of 
~0.30 magnitudes in J and of ~0.10 magnitudes in the H and Kj 
bands for the NACO-to-NACO registration, while the scatter in- 
creases to ~0.20 magnitudes for the Gemini-to-NACO registra- 
tion. Hence, we take these values as representative of the random 
errors of our relative photometry. By adding them to the formal 
errors we find that only the H band variability (1.07+0.26) can 
be considered formally significant, and it is consistent with that 
measured in the J band (0.97+0.43) and, marginally, also with 
the Kj band one (0.40+0.16). 

3.2. 1RXSJ1 708-4009 

Fig. [3] shows the NACO A'j-band image of the field around the 
Chandra position of IRXS J1708-4009 (RA 11^ 08™ 46^9, DEC 
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Fig. 4. IRXS J1708-4009: Colour-magnitude diagrams for the objects in the field of view. Sources detected within or close to the 
X-ray error circle (see Fig. 13) are shown in red and labeled. 



-40° 08' 52'.'4). The size of the X-ray error circle has been esti- 
mated considering the un certainty on the Chandra position, that 
is 0'.'7 (llsrael et al.ll2003 l). combined with the overall accuracy 
of our astrometric caUbration (see §2.2), yielding a final uncer- 
tainty of r.'l at 99% confidence level. The magnitudes of all 
the objects detected within the X-ray error circle are reported in 
Tab. [3] No other object is detected down to the 3cr limiting mag- 
nitudes of 22.98, 22.20, 22.26 in the J, H, Kj bands, respectively. 
In the L' band no object is detected within the X-ray error circle 
down to 17.8 (the deepest limit ever obtained for an AXPs in 
this band). Fig.g] shows the (K,, J-K,) and (H, H-K,) CMDs of 
all the objects detected in the field. The two CMDs show a fairly 
scattered sequence. Despite previous claims (Israel et al. 2003; 
Safi-Harb & West 2005; Durant & van Kerkwijk 2006), none of 
these objects show significant (>3cr) IR variability to justify a 
safe identification as the IR counterpart of IRXS J1708-4009 
(see also ®. 

3.3. 1E 1841-045 

An updated determ ination of the IE 1841-045 position 
dWachter et al.ll2004l) was obtained with Chandra and gives RA 
Igh 19s 336^ DEC -04° 56' 10'.'83, with a 3cr uncertainty 
radius of 0'.'9. However, the latter uncertainty is based only on 
five 2MASS reference stars used for the boresight correction 
(Wachter et al. 2004). We decided to use a conservative nom- 
inal Chandra positional uncertainty of 0'.'7 at 90% confidence 
level. Fig.|5] shows a revised finding chart of the field around 
the updated X-ray position. The size of the X-ray error circle, 
r.'2 (99% c.l.) accounts for the overall accuracy of our astromet- 
ric calibration (see §2.2). The magnitudes of all objects detected 
within or close to the X-ray error circle are reported in Tab.|4] No 
other object has been detected down to the 3cr limiting magni- 
tudes of 21.54 and 21.00 for H and K', respectively. All objects 
have been detected in both the H and K' bands, except for #9 and 
#10. Fig. |5] shows the position of these objects in the (K',H-K') 
colour-magnitude diagram compared to other objects detecte d 
in the field. We note that object #19 of iMereghetti et all (l200lh . 
originally considered a potential candidate on the basis of the 



ROS AT X-ray position, lies more than two error radii away from 
the updated Chandra position. The CMD shows a well-defined 
main sequence extending up to K' ~ 12, with the candidates oc- 
cupying the fainter region and stretching over a broad range of 
colours. 

By comparing our photometric results with those reported in 
Uterature for IE 1841-045 (Durant 2005), it is evident that our 
source #9 (source B in Durant (2005)), shows a >4cr variabil- 
ity (AKs =1.30±0.24). As a comparison, the average AKs and 
AH between our photometry and that of Durant (2005) for 9 ob- 
jects in the field of view is ~0.2 magnitudes. We then tentatively 
propose this object as the IR counterpart to this AXR 

3.4. SGR 1900+14 

Fig. IS] shows the NACO band image around the 

SGR 1900+14 radio position (RA 19'^ 07™ 14^33, DEC 
-1-09° 19' 20'.'!, with a 1 cr uncertainty of 0'.'15 in each coor- 
dinate; Frail et al. 1999). Our astrometry yields an overall 3cr 
position uncertainty of 0'.'81, which comp ares well with the on e 
of 0'.'79 (99% confidence level) quoted bv lKaplan etal] (|2002|) . 
We note that the PSFs of field objects are slightly elongated in 
the NACO images, and one of the objects detected within the 
radio error circle (our #3 and object C of Kaplan et al. (2002)) 
looks apparently extended. However, its PSF is consistent with 
that measured for most field objects. In this case, photometric 
calibration was performed using as reference the magnit udes of 
a few objects det ected around the target position (iKaplan et alJ 
I2OO2I iFrail et al.l [1999 ) (see Fig.©. The three new objects 
detected in our NACO images are also marked (# 7, 8 and 9). 
In particular, object #7 falls very close to the position of #3, 
which prompted us to use PSF photometry rather than aperture 
photometry as done, instead, for the other NACO observations. 
The magnitudes of all the detected objects are listed in Tab. 
|5]No other object is detected close or within the radio position 
error circle down to a 3cr limiting magnitude of = 22.60 
and 22.25 for the March 31 and for the July 20 observations, 
respectively. We note that among our newly detected objects, #7 
(K5 ~ 19.7) is the only new source within the 99% error circle. 
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Fig.5. Left panel: IE 1841-045 : CFHT K'-band image of fie ld, with the ir2 radius (99% confidence level) X-ray error circle 
overlaied. Objects detected at > 3cr are labeled. Object #19 of iMereghetti et alj (1200 1 [) is labeled as a reference. Right panel: K', 
H-K' CMD for IE 1841-045 field. Objects detected within or close to the X-ray error circle (left) are shown in red and labeled. 



and it displays some evidence of variability (A K, =0.47+0.11) 
between the two epochs. Unfortunately, its proximity to object 
#3 prevents a better determination of the actual magnitude 
variation, even through PSF photometry. As we did in §3.1, 
we evaluated the effects of random errors on our relative 
photometry by measuring the magnitude scatter of a number 
of test stars of brightness comparable to object #7 and selected 
close to the field center in order not to introduce biases in our 
PSF photometry due to the position-dependent NACO PSF. We 
found that for the test stars the scatter is « 0.1 magnitudes. Even 
adding the random error we found that the measured variation 
is still significant, although only at the ~ 3cr level. Based on 
variability, we then propose object #7 as a tentative IR candidate 
counterpart to SGR 1900-H14 . 

4. Discussion 

Before discussing our results, we would like to clarify here how 
we define a reliable identification of an IR counterpart to an AXP 
or SGR. All of the current well defined counterparts have been 
established either through a very accurate positional coincidence 
(no other candidates in the positional error circle) or by the de- 
tection of significant IR variability. In some cases, what are usu- 
ally defined (not very clearly) as "strange" colours with respect 
to the other stars in the field of view has been used as a method 
to claim the detection of a counterpart. However, we prefer not 
to adopt this "colour" criterion because we find it rather mislead- 
ing. In fact, i) almost all counterparts reliably identified through 
IR variability detection or positional coincidence have colours 
not much at variance with some of the their field stars, and ii) 
our current lack of knowledge of the the exact physics behind the 
IR emission of AXPs and SGRs does not allow any prediction 
on their IR spectra. Another criterion, often used to propose or 
strengthen an identification is the X/IR flux ratio. We believe this 
method rather rough and not constraining at all for the following 
reasons: i) AXPs and SGRs are variable both in the X-ray and 
IR bands, and we do not have yet a clear idea on the connection 
between these two bands, ii) there are no theoretical predictions 



with respect to the X/IR ratio that an AXP or SGR is expected 
to have, and iii) with the exception of their IR variability, the IR 
counterpart to these neutron stars are, at our current knowledge, 
consistent with the IR emission of very low mass stars, which if 
laying in the source eiTor circle will then result in the same X/IR 
ratio. This is why as long as significant variability is not detected, 
if more than one faint IR object is present in the positional circle, 
any proposed IR conterpart needs a further confirmation. 

4.1. XTEJ1810-197 IR variability 

IR observations with the VLT-NACO camera of XTEJ1810- 
197 performed after th e outburst, unveiled its IR counter- 
part jlsrael et al.l l2004bl) . Follow-up observations performed 
six months later detected IR variability from the candidate 
counterpart (jReaetal .12004 . with a flux decrease by about a fac- 
tor of 2. Simultaneous X-ray observations over the same period 
revealed a similar X-ray flux decrease between the two epochs 
(Rea et al. 2004; Gotthelf et al. 2004). This suggested that the 
coiTelated IR/X-ray variability is a characteristic of this source, 
as expected if the IR emission were dominated by the reprocess- 
ing of the X-ray emission from a fossil disk around the magnetar 
(Perna et al. 2000; Perna & Hernquist 2000). However, our re- 
cent IR observations show that such correlation is not as obvious 
as thought before (see also Camilo et al. 2007b). While the IR 
flux is indeed variable, it does not follow the same monotonic 
decrease of the X-ray flux observed in the post-outburst phase 
(see FiglHand Gotthelf & Halpern 2007). Instead, the IR vari- 
ability appears more similar to the one observed in radio band 
(Camilo et al. 2006, 2007c). Unfortunately, the IR observations 
are too sparse for any firm conclusion to be drawn regarding 
a possible connection between the radio and the IR variabil- 
ity. Furthermore, we note that the onset of the radio emission 
(January 2005; Camilo et al. 2006) occurred after our IR obser- 
vations. 

In the fossil disk scenario, the XTEJ1810-197 IR variabil- 
ity can be hardly explained if the IR emission were dominated 
by reprocessing of the X-ray radiation from the pulsar. In this 



8 



Testa et al.: Adaptive optics near-IR observations of magnetars 



Id. 


RA (°) 


DEC C ) 


H 


K' 


1 


280.330811 


-4.936699 


20.70± 0.31 


19.53± 0.15 


2 


280.330444 


-4.936720 


20.62± 0.29 


20.95± 0.49 


3 


280.330597 


-4.936695 


19.58±0.10 


18.19± 0.03 


4 


280.330658 


-4.936555 


18.48±0.03 


18.56± 0.05 


5 


280.330475 


-4.936538 


19.01±0.06 


19.86± 0.18 


6 


280.330627 


-4.936442 


19.42± 0.08 


19.24± 0.11 


7 


280.330933 


-4.936414 


17.88± 0.02 


I7.07± 0.02 


8 


280.330353 


-4.936280 


20.69± 0.29 


18.31±0.04 


9 


280.330760 


-4.936186 


>21.54 


20.47 ± 0.28 


10 


280.330670 


-4.936334 


>21.54 


19.93 ±0.16 



Table 4. IE 1841-045 : H and K' magnitudes for all the objects 
detected within 17 1 of the source X-ray position. Magnitude er- 
rors are given at Icr confidence level. 



Id. 


RA(°) 


DEC (°) 


K(2006.03.31) 


K(2006.07.20) 


1 


286.81000 


9.32198 


19.36±0.09 


19.39±0.06 


2 


286.80949 


9.32190 


18.14±0.07 


18.05±0.05 


3 


286.80957 


9.32212 


17.31±0.07 


17.25±0.05 


4 


286.80948 


9.32216 


18.42±0.07 


18.48±0.05 


5 


286.80989 


9.32189 


20.10±0.11 


19.96±0.08 


6 


286.80974 


9.32220 


20.73±0.26 


20.85±0.22 


7 


286.80963 


9.32215 


1 9.2 1 ±0.08 


19.68±0.08 


8 


286.80997 


9.32206 


20.86±0.19 


20.55±0.13 


9 


286.81016 


9.32210 


20.59±0.20 


20.46±0.15 



Table 5. SGR 1900+14: K magnitudes of the March and July 
2006 observations, for all objects within 1'.'6 from the position 
(see text for details). Errors on the magnitudes are given at Icr 
confidence level. 



case in fact the IR flux should decrease as the X-ray flux drops 
(see Rea et al. 2004 for details), and no increase in any of the 
IR bands would be expected. If, on the other hand, the IR flux 
were dominated by the disk emissivity resulting from viscous 
dissipation (which could be the case if the disc inner radius is 
farther out), then flux variability would be the result of a varia- 
tion in the disk mass inflow rate. Given the non-correlated IR and 
X-ray variability we observe, this latter scenario could be more 
viable. However, while a variability in the IR emission is easily 
predictable as resulting from a mass inflow rate variation, this 
should take place similarly in all the IR bands. As noted in ^3.11 
we see a hint for an IR-band dependent variability, although no 
significant conclusion can be drawn on this aspect yet. 

Within the magnetars scenario, a few different models have 
been proposed for the production of the IR radiation (e.g. Heyl 
& Hernquist 2005; Beloborodov & Thompson 2007), but none 
of them makes (yet) specific predictions for the IR variability 
that would permit a direct comparison with our data. 

IR variability has been observed for other magnetars, 
i.e. IE 1048.1-5937 (Israel et al. 2002; Tam et al. 2007), 
4U0142-H614(Hulleman et al. 2004), SGR 1806-20 (Israel et 
al. 2005), and 1E2259-H586 (Tam et al. 2004). Unfortunately, 
with the current data no firm conclusion can be drawn yet about 
the IR variability in connection with variabilities in other observ- 
ing bands, except for 1E2259+586 and SGR 1806-20 (Tam et 
al. 2004; Israel et al 2005), where a correlation with the source 
bursting activity and X-ray flux enhancement has been detected. 

4.2. Other sources: 1E 1841-045, SGR1900+14 and 
1RXS J 1708-4009 

No IR counterpart had been claimed so far for IE 1841-045 and 
SGR 1900+14, likely because of the crowding of the field in 
which these magnetars happens to lay. By comparing our photo- 
metric results with those reported in literature for IE 1841-045 , 
it is evident that there is one object, our source #9 (source B 
in Durant (2005)), with magnitudes similar to the other AXP 
counterparts, and showing a >4cr variability. Furthermore, com- 
paring our two observations of SGR 1900+14 we found a ~3cr 
variability in source #7 (not detected by Kaplan et al. (2002) 
because it was too faint for their limiting magnitudes). We there- 
fore consider these objects as possible IR counterparts to these 
two neutron stars. 

The new deep observation of IRXS J 1708^009 show 20 
sources within and in the vicinity of the 99% X-ray error circle 
(see Fig. [3). None of these objects showed significant variability 
with respect to other observations in the literature, which could 
help identifying them as reliable magnetar counterparts. 



Note that the magnitudes of both the previous candidates 
(our #1 and #3; Israel et al. 2003, Safi-Harb & West 2005, Durant 
& van Kerkwijk 2006) would make IRXS J 1708-4009 much 
brighter in the IR band than any other AXP (see e.g. Israel et al. 
2004a). Furthermore, all the other fainter candidates we report 
here, laying within the source positional errors (see Fig. 3), have 
brightness in better agreement with th e IR emission properties 
of AXPs than both objects #1 and #3. iDurant & van Kerkwiilj 
(2006) proposed source #3 as a possible candidate on the basis 
of IR variability, the strange colours and the resulting X/IR ratio 
of this object. We believe that given the large number of faint 
IR sources laying within the IRXS 11708^009 positional er- 
ror circle, and the rather bright magnitudes of source #3, neither 
the ~ 2.5cr IR variability, nor its colours (see Fig.lUi or its X/IR 
ratio can reliably tight this object to the AXP (see also above). 
We would like to stress that we can not exclude any of these 
candidates, but that a firm identification of an IR counterpart to 
IRXS J 1708-4009 is far from being a settled issue. 

This AXP has been recently observed to have a vari- 
able X-ray emission (Rea et al. 2005; 2007b; Campana et 
al. 2007), which, if correlated somehow with the IR as for 
IE 2259+586 (Tam et al. 2004), would imply a variable IR 
counterpart. In particular, the X-ray variability observed for 
IRXS J 1708-4009 in X-ray observations sparse over several 
years was of the order of 50%, which would roughly imply (if 
X-ray and IR are directly correlated) an IR variability of ~ 0.5 
magnitudes, not detected from any of the objects near the posi- 
tion of IRXS 11708^009. However, the lack of simultaneous 
IR and X-ray observations of this source prevent us from draw- 
ing any firm conclusions about possible correlated variability. 

5. Conclusions 

The IR counterpart to XTEJ1810-197 is confirmed to be vari- 
able in time (as previosuly proposed by Rea et al. (2004)). 
However, despite previous claims, its variability might be more 
similar to its radio behavior than to its X-ray variability, which 
appears to decay smoothly after the outburst. Simultaneous IR 
and radio observations of this AXP are needed to reliably assess 
this possibility. 

To date, IR counterparts to AXPs and SGRs have been 
confirmed only for XTEJ1810-197 (Israel et al. 2004b), 
IE 1048. 1-5937 (Wang & Chakrabarty 2002; Israel et al. 2002), 
1E2259+586 (Hulleman et al. 2001), 4U0142+614(Hulleman 
et al. 2000) and SGR 1806-20 (Israel et al. 2005; Kosugi et al. 
2005). For the AXP IE 1841-045 and SGR 1900+14, we pro- 
pose here two possible candidates based on the detection of IR 
variability. For all the remaining AXPs and SGRs we still miss 
a candidate or confirmed IR counterpart (see e.g. Wachter et al. 
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Fig, 6. SGR 1900+14. Zoom-in of the area of the Kj band im- 
age of 2006 March 31st around the radio position of Frail et 
al. (1999). The error circle (radius O'.'Sl at 99% confidence 
level) is shown. Two large ghosts are visible atop the position 
of SGR 1900+14 , hiding star E of Kaplan et al. (2002). 

2004; Durant & van Kerkwijk 2005, 2007; Muno et al. 2006; 
Gelfand & Gaensler 2007). 

Ours and others results on optical and IR observations of 
AXPs and SGRs, show that up to now we are far from having 
an overall picture of the optical and IR behavior of these neutron 
stars. Further observations, possibly simultaneously on a wide 
energy range, are needed to refine current theoretical models and 
shed light on AXPs' and SGRs' optical and IR emission mecha- 
nisms. 
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